Electronic Properties of Armchair Carbon Nanotubes : Bosonization Approach 
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The phase Hamihonian of armchair carbon nanotubes at half-filling and away from it is derived 
from the microscopic lattice model by taking the long range Coulomb interaction into account. We 
investigate the low energy properties of the system using the renormalization group method. At 
half-filling, the ground state is a Mott insulator with spin gap, in which the bound states of electrons 
at different atomic sublattices are formed. The difference from the recent results [Phys. Rev. Lett. 
79, 5082 (1997)] away half-filling is clarified. 
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Single wall carbon nanotubes (SWNTs) with diame- 
ters of a few atomic distances and lengths of several mi- 
crometers can be considered as the ultimate miniatur- 
ization of metallic wires Recent experiments have 
demonstrated electron transport through individual ||] 
and multiple SWNTs as well as provide evidence of 
strong Coulomb interaction in these systems. The one- 
dimensional nature of the low-energy electronic states in 
the nanotubes together with the interaction of electrons 
should result in a variety of correlation effects due to the 
non- Fermi liquid ground state of the system . 

Very recent transport spectroscopy data by Tans et. al. 
1^ on spin polarization of an individual SWNT can not be 
explained by the constant interaction model and suggests 
the interpretation in terms of electron correlations. This 
result, however, was not confirmed by experiments on 
ropes of SWNTs which fit the constant interaction 
model remarkably well. 

Experimental progress urges the development of a the- 
ory of electron correlations in SWNTs. For a model on- 
site [0 and on-site plus nearest neighbor Q interactions, 
metallic armchair SWNTs become Mott insulator at half- 
filling, whereas upon doping they exhibit superconduct- 
ing fluctuations. The realistic long-range Coulomb in- 
teraction was considered in Refs. Kane, Balents 
and Fisher discussed the effects of the Coulomb 
interaction in finite size armchair nanotubes (ANs) in 
terms of the Tomonaga-Luttinger low energy theory. The 
most important part of the forward scattering was incor- 
porated into the Tomonaga-Luttinger-like Hamiltonian, 
whereas the other types of scattering were treated as per- 
turbations. Egger and Gogolin (EG) developed an effec- 
tive low energy theory of ANs starting from a microscopic 
model ||] , which accounts for all types of scattering pro- 
cesses. They derived a bosonic phase Hamiltonian and 
discussed possible ground state away from half-filling. 

In this Letter we derive the phase Hamiltonian of ANs 
and evaluate its parameters from the microscopic lat- 
tice model. The difference between our Hamiltonian and 
that by EG stems from the distinction in the form of ki- 
netic term and the use of oversimplified approximation 



for 2fci^-component of scattering amplitudes in Ref. 0, 
which overlooks C3 symmetry of the microscopic lattice 
model. The renormalization group (RG) method is ap- 
plied to the Hamiltonian, and the low energy states are 
investigated. At half-filling the ground state is the Mott 
insulator with spin gap, in agreement with the conclusion 
of Hubbard- like models In this state, the electronic 

bound states are formed between the different sublat- 
tices. Away from half-filling we predict gaps for both 
symmetric and antisymmetric spin modes, in contrast to 
the result of Ref. (|] for the case of equal amplitudes of 
intrasublattice and intersublattice forward scattering. 

We start from the standard tight-binding single parti- 
cle Hamiltonian pd] | for pz electrons on the honeycomb 
lattice (inset of Fig. 1), 



(1) 



Here ap^s{k) are the Fermi operators for electrons at the 
sublattice p — ± with the spin s = ± and the wavevector 
k — {kx.,ky). The matrix elements are given by £^{k) = 
_^(-g-ife„a/\/3_|_2eifeH<i/2v^(,Qg^^Q,^2), t being the hopping 
amplitude between neighboring atoms. The eigenvalues 
of the Hamiltonian vanish at two points of the Brillouin 
zone, k — {aKo, 0) with a = ± and Kq = 47r/3a, which 
constitute the Fermi surface of a graphite layer |11 ^. 

We consider the armchair {N,N) SWNT parallel to 
the X axis so that the wrapping vector w = N(d-\. + 
a_) points in y direction (inset of Fig.l). In this case 
the Fermi points lie on the allowed quantized transverse 
wavevector fc^ = for any N. Expanding Eq. (^ near 
the Fermi points to the lowest order in q = k^ — aKo 
and introducing slowly varying Fermi fields ^^pasi^) = 



Eg e'l^'ap.siq + aKo, 0), we obtain 



(2) 



p,a,s 



= ^/3ta/2 w 8 X 10^ m/s being the Fermi velocity. It 
should be noted p| that Eq. (|) differs from Eq. (2) of 
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Ref. . Despite identical energy spectra of both Hamil- 
tonians, the phase relations between the components of 
an eigenfunction at the two sublattices are different. The 
consequences of this fact will be discussed later on. 

Following EG, the interaction term of the Hamiltonian 
reads, 

^»*-^(i)'EE E (3) 

1,1' pp' ai...a4 ss' 

with xi — la/2. The effective ID interac- 
tion between the sublattices p and p' , Upp'{xi), is 
the average of the Coulomb potential U{x,y) — 

/{k^ flg + + AB? s\v?{y /2R)} over the nodes of a 
sublattice along the y direction, 

Upp' {xi) = -j^^U (^xi,a\/2,{n + A. 



■pp' ) j J 



(4) 



with Appi = mod(Z, 2)/2 + Sppi /3. Here k is an effective 
dielectric constant of the system (the estimate Q for the 
parameters of the experiment |^ gives k — 1.4) and oq ~ 
a characterizes the radius of pz orbital wave function in 
the graphite plane. 

The Hamiltonian Timt can be separated into the "for- 
ward scattering" Ho + Ti-f (ai = a^, 012 = a^) and 
" backscattering" Tib + Ti-t' (ai ~ —ce2 0^3 = —cti) 



Ho — 
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5V(Q) 



dxp^{x)^ 



(5) 



Hh 



Vpp{2Ko) 



^ j da;i/;t^3Vlpa's'V'-pa's'V'pas, (6) 

\a.' ss' 

E / da^V'pasV'p'-as'V'p'as'V'p-a*, (7) 



pp' ass 



Hb' — — 



SV{2Ko) 



pass' 



where p = J2pas i^lasi^pas is the total electron density, 
SV{0) = VppiO) - Vp_p{0), and 6Vi2Ko) = Vpp{2Ko) - 
Vp_p{2Ko), with Vpp,{q) = (a/2)^,e^«-'C/pp,(a;,) being 
the Fourier transformed interaction. 

The forward scattering (||) has the strongest ampli- 
tude, Vpp{0) = 2e'^/Kln{Rs/R), where Rs ~ min(L,D) 
characterizes the screening of the Coulomb interaction 
due to a finite length L of the AN and / or the presence 
of metallic electrodes at a distance D . The ampli- 
tude Vpp{0) is relatively insensitive to the choice of Rg 
due to logarithmic dependence. From Eq. (^), one sees 
that the amphtudes 6V{0) and Vpp{2Ko) decay as 
for R ^ a. It should be noted that the matrix element 
NVp-p{2Ko) vanishes identically in case of a graphite 
plane {R ^ oo) due to the C3 symmetry of the lat- 
tice. For this reason, Vp-p{2Ko) is much smaller than 



5V{Q) and Vpp{2Ko). All the matrix elements decrease 
with increasing oq. Numerical evaluation for ag = a/2 



and i? > a gives dV{0) = 0.21, Vpp{2Ko) 



0.60, 



Vp^p{2Ko) = 9.4x10"'' in units ofae'^/2nKR {Vp-p{2Ko) 
is estimated for TV — 10). This result shows that the ap- 
proximation Vpp{2Ko) = Vp-p{2Ko) used in Ref. ^ is 
questionable. 

In order to bosonize the Hamiltonian H ~ Hk + Hint, 
Eqs. (^), (||)-(||), we diagonahze the kinetic term Eq. (||) 
by the unitary transformation 



Aas = {i>+as + arip^as)/ %/2, 



(9) 



which maps the basis of atomic sublattices {p — ±) to 
the basis of right- and left-movers (r = ±). This trans- 
formation is different from one by EG due to the different 
form of the kinetic term. 

We bosonize the Fermi fields V'ras, 



, ^r,a.s 
Ipras = ,—— exp 



WqpX + — {Oas + r(t)as} 



(10) 



and decompose the phase variables Oas,4'as into sym- 
metric S = + and antisymmetric 6 = — modes of the 
charge p and spin a excitations, 9as = (^p+ + s9cr+ + 
aOp-. + asOrj^ and 4>as = 4>p+ + s4>a+ + cx-i'p- + as^cr-- 
The bosonic fields satisfy the commutation relation, 
[Ojs{x),(j)j'S'{x')] = i{-K /2)sign{x - x')8jj' 5ss' . The Majo- 
rana fermions rjras are introduced to ensure correct an- 
ticommutation rules for different species r, a, s of elec- 
trons, and satisfy [riras,ilr'a' s']+ = 25rr'5aa'5ss' ■ The 
spin-conserving products rh-asVr'a's in the Hamiltonian 
H can be represented as A++{r, a, s) = rjrasVras = 
1, ^+_(r, a, s) = rjrasVr-as = iacTj,, A_+(r, a, s) = 

TJrasTl-ras = ^raa^,, and A (r, a, s) = rirasV-r^as = 

—vray with the standard Pauli matrices Ui (i = x,y,z). 
The quantity qp — T^n/A is related to the deviation n of 
the average electron density from half-filling, and can be 
controlled by the gate voltage. 

The bosonized Hamiltonian has the form, 

^ = E EI7 Ux[K-l{dJ,sf+K,s{d..cp,sf} 



2(TTa)'^ J 

[SV{0) - 2V^(2i^o)] cosiAqpx + 29 p+) cos2e. 
-5V{0) cos{AqFX + 26 p+) cos 2dp- 
+5V{0) cos{4qFX + 29 p+) cos 20^- 
~ [5V{Q) ~ 5V{2Ko)] cos 29 p_ cos 29^,- 
+5V{Q) cos 26*^+ cos26i^_ 
-5V{Q) cos 26*^+ cos26ip_ 
-2V{2Ka) cos{4qFX + 29 p+) cos 2(/)^„ 
+2V{2Ko) cos 261^+ cos 20<,_ 
+5V{2Ko) cos 29 cos 20„_ 
+5V{2Ko) cos 29„^ cos 2<l)„^}, 



(11) 
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Vjs = vo^/ AjsBjs and Kjs = y^Bjs/Ajs being the veloc- 
ities of excitations and exponents for the modes j, S. The 
parameters Ajs, Bjs are given by 



A„. = 1 



Bp 
A, 

Ap 



Bn 



4F(0) 6V{Q) V{2Ko) SV{2Ko) 

ttvq Attvq 2t:vo -iTTVQ 

_ ^ ^ dVjO) ^ Vi2K„) 6Vi2Ko) 

^ Attvo 27ruo Anvo ' 

5V{0) V{2Ko) 6V{2Ko) 

Airvo 2t:vo Attvq ' 

SV{0) V{2Ko) SV{2Ko) 



= A„.^l- 



= B„^ = l 



5V{0) _ V{2Kq) 
47rwo 2'kvq 



SV{2Ko) 
Attvo 



(12) 
(13) 
(14) 
(15) 
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with V{q) = [Vpp{q) + Vp-p{q)]/2. The sublattice- 
independent forward scattering V{0) strongly renormal- 
izes the exponent for the symmetric charge mode, Kp^ « 
0.2 [0, whereas for the other modes the interaction is 
weak, Kjs = 1 + 0{a/R) §. 

We now compare the Hamiltonian (|ll]) with that de- 
rived by EG Q away from half-filling. In this case, the 
non-linear terms containing the misfit parameter qp can 
be neglected due to the breakdown of the momentum 
conservation. Despite the equal forward scattering parts 
of both Hamiltonians, there is difference in the backscat- 
tering parts. Namely, the Hamiltonian by EG can be 
obtained from ours by substituting V{2Ko) and 
6V{2Ko) 2V{2Ko). The absence of cos 26*^+ cos 20<,_ 
term in the Hamiltonian by EG originates from the use 
of the approximation 5V{2Kf)) = 0, whereas the differ- 
ence in the coefficient in front of the three other non- 
linear backscattering terms stems from the difference in 
the unitary transformation (P). 

The low energy properties of the Hamiltonian Eq. 
( |rT| ) can be investigated by the renormalization group 
(RG) method. The RG equations can be derived by 
assuming scale invariance of the correlation functions 
(exp {i (^,5(xi, Ti) - e,s{x2,T2))}) 11. At half-filling, 
qp = 0, we obtain 



{Kp+y 

{Kp-y 
(K.-y 



-{Kl 

(l/8)(2/? 
{2 - {Kp+ 4 

(2/2^6 + 2/32/5 

{2 - {Kp+ 4 
(2/12/6 + 2/32/4 
{2 - {Kp+ 4 



2/1 



yl 

-vl 
-vl- 
-vl 
vl) 



2/3^ 

- 2/6 ^ 
-vl) 



-2/?) 

-vl) 
vl) 



(2/i)' 
(2/2)' 
(2/3)' 

(2/4)' = {2-(i^p-+i^.-)}2/4 



AW)} 2/1 
-l-2/72/8)/4 

Kp-y\vi 

-l-2/72/9)/4 

Ka-)}VZ 



(2/12/5 +2/22/4)74 



(17) 
(18) 
(19) 

(20) 
(21) 
(22) 
(23) 



- (2/22/3 + 2/52/6)/4 , 
(2/5)'-{2-(i^.+ +if.-)}2/5 

- (2/12/3 + 2/4'2/6)/4 , 
(2/6)'-{2-(i^.+ +i^p-)}2/6 

- (2/12/2 + 2/42/5 + 2/82/9)/4 , 
{y^y = {2-{Kp^^\|K,J)\y^ 

- (2/12/8 + 2/22/9)74 , 

(2/8)' = {2-(X.+ + l/if.-)}2/8 

- (2/12/7 + 2/62/9)/4 , 

(2/9)' = {2-(Xp_ + l/X,_)}y9 

- (2/22/7 + y&y%)lA , 



(24) 

(25) 
(26) 
(27) 
(28) 
(29) 



where ()' denotes d/d^ with d£ = dln(a/a) (a is the new 
lattice constant). The initial conditions for the Eqs. (|l7|)- 
(H) are K,s{^) = K,^, y^ = [5V{Qi) - 2V{2Ko)]/ {ttv^), 
2/2 = -2/3 = -2/5 = 2/6 = -SV{Q)/{t:v^2), y^ = -[6V{Q) - 
5V{2Ko)]/{^vo), 2/7 - -2/8 - -2V{2Ko)/{^vo), and 
j/g = 5V {2Kq) / {-KVa) . In deriving the RG equations, 
the non-linear term cos 20o— cos 2(/)o— is omitted because 
this operator stays exactly marginal in all orders and is 
thus decoupled from the problem [|l5|. The RG equa- 
tions away from half- filling can be obtained from Eqs. 
(|r7|)-(p9|) by putting yi, y2, 2/3 and yj to zero. Hereafter 
we concentrate on the case TV = 10, k = 1.4, = 100 
nm and oq = a/2 where the initial values of the RG 
parameters correspond to the estimates given below Eq. 



Away 
quantities 



from 

K„^ and 



half-filling, the 
i rcnormalize to zero and 



cos 26*^+ cos 2(/)cr-and 
00). As a result, the 



the coefficient of cos 26a-+ cos 26'p_ 
cos 20p_ cos 20cr-) tends to —00 (00). As a result, tl 
phases 6*^+, 6'p_ and (pa-^ are locked at (6*0-+, ^'p-, 4>a--) 
(0, 0, 7r/2) or (7r/2, 7r/2, 0) so that the modes a± and p— 
are gapped. In this case, the asymptotic behavior of the 
correlation functions at x — *■ 00 is determined by the cor- 
relations of the gapless p+ mode, ^e'"^''+^^^e~™^f+(°)) ~ 



(n ■ 



1 and 2 correspond to 2qp and Aqp density waves and 
TO = 1 for a superconducting state). Since -ftrp+ ~ 0.2, the 
2qp density wave correlations seem to be dominant. How- 
ever, we found that the correlation functions of any 2qp 
density wave decay exponentially at large distances due 
to the gapped modes. We therefore are looking for the 
four-particle correlations. The Aqp density waves domi- 
nate over the superconductivity for Kp+ < 1/2 |]l6| , |r7| . 
Such density wave states are given by the product of the 
charge n± (x) or spin S± (x) densities at different sublat- 
tices. 



S+{x)S- (x) 



■ cos{Aqpx + 29p+) 



2(7ra)2 

X (2 cos 29^+ + cos 20<^_ - cos 26'p_), (30) 

— i-^ cos{Aqpx + 20 p+) 

X (2 cos 29^+ - cos 20<,_ + cos 26'p_), (31) 
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where we neglected the unlocked phases 0o-_|_ , 0p_ , 
whose contribution decays exponentially at large dis- 
tances. Substituting the values of the locked phases 
we observe that vanishes, and the dominant 

state is the Aqp spin density wave with correlation func- 
tion {S+{x)S-{x)S+iO)S-{0)) ^cos4qFx/x^^''+. 

The modes ct± and p— remain gapped also in the limit 
6V{0) = 0. In this case EG have obtained that the sym- 
metric modes, p+ and a+, are gapless, whereas the p— 
mode is gapped and a— mode separates into the gapless 
and gapped parts. The result by EG follows from the 
special dual symmetry 9a-- 4>a- of the Hamiltonian 
and the absence of non-linear terms in the p+ and a+ 
sectors. Both these factors are lacking in Eq. (^. On 
the other hand, the result by EG for a finite value of 
6V{0) is qualitatively the same as ours. 

At half- filling the solution (Fig.l) of the RG equations 
(p7|)-(p9|) indicates that the phase variables 6'p+, 0^+, Op-, 
and ^CT- arc locked and the all kinds of excitation are 
gapped. In other words, the ground state of the half- 
filled AN is a Mott insulator with spin gap. The same 
conclusion has been drawn from the model with short 
range interactions (7|J|]. The locked phases are given by 
{9p+,9a+,9p-,(pa-) = (0,0,0,0) or (^2, 7r/2, ^/2, 7r/2) 
since the coefficients tend to — oo for the first, second 
and 6-9-th non-linear terms in Eq. (|Tl|). The averages 
{n+{x)n-{x)) and {S+{x)S-{x)) are both finite, which 
indicates the formation of bound states of electrons at 
different sublattices. 

The observability of the Mott insulating behavior of 
ANs depends critically on the magnitude of the gap Ap+ 
in the p+ mode, which is estimated by the self-consistent 
harmonic approximation as , 

Ap+/(2vp+a-^) 



5V{Qf 



V{2KoY - dV{0)V{2Ko) 



(32) 



Using the value of the matrix elements calculated nu- 
merically, the charge gap for iV = 10 is estimated as 
^ lOOK for ao = a/2 and ^ lOK for ao = a (in case of 
ao = a, SV{0) = 5.6 x lO'^ and V{2Ko) = 6.9 x lO'^ 
in units of ae^ /2ttkR). The resistivity of ANs shows 
power law temperature dependence, p ~ T^^p+~^/A^^ at 
high temperatures T 3> Ap+ and increases exponentially, 
p cx exp(Ap+/T), at T < Ap+ |l0|. The temperature de- 
pendence of the resistivity at half- filling is a characteristic 
signature of the Mott transition. We conjecture that this 
signature can be best detected in multiprobe transport 
measurements 119(1. 
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FIG. 1. Solutions of the renormalization group equations 
for Kp+ , Kcr+ , Kp- and K^- . The initial values of parameters 
correspond to the estimates given in the text below Eq.(8) 
for A = 10, K = 1.4, Rb — lOOnm and ao — a/2. Insert: 
The honeycomb lattice of carbon atoms. Here a± are the 
two primitive Bravais lattice vectors, |a±| = a. The hexagon 
shown by thick line is the unit cell and the black (white) circle 
denotes the point at p — +( — ) sublattice. The x axis points 
along AN. 
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